Escape of a Uniform Random Walk from an Interval 



in 

O 

o 

O 

Q 

(N 
(N 



S3 

c3 

c$ 
+-> • 
ctf ■ 

O ■ 

• i-H . 

C/3 . 
>v 
43 ■ 



<n : 

(N . 
(N ■ 

in : 
o . 

o : 
>>: 

Oh' 



X 



T. Antal 1 and S. Redner^Q 

1 Center of Polymer Studies and Department of Physics, 
Boston University, Boston, Massachusetts, 02215 USA 

We study the first-passage properties of a random walk in the unit interval in which the length 
of a single step is uniformly distributed over the finite range [—a, a]. For a of the order of one, the 
exit probabilities to each edge of the interval and the exit time from the interval exhibit anomalous 
properties stemming from the change in the minimum number of steps to escape the interval as a 
function of the starting point. As a decreases, first-passage properties approach those of continuum 
diffusion, but non-diffusive effects remain because of residual discreteness effects. 

PACS numbers: 02.50.C2, 05.40.Fb 



I. INTRODUCTION 

Consider a discrete-time random walk in which the 
length of each step is uniformly distributed in the range 
[—a, a]. We term this process the uniform random walk 
(URW). The walker is initially at an arbitrary point x in 
the unit interval [0,1] and the endpoints are absorbing. 
For the URW, an absorbing boundary is defined such that 
if the walk crosses an endpoint of the interval, the walk 
is trapped exactly at this endpoint. We are interested in 
the first-passage properties of this URW. 

One motivation for this study comes from the prob- 
lem of DNA sequence recognition by a mobile protein 
P, ■ The protein molecule typically seeks its target by 
a combination of diffusion along the DNA chain, and also 
detachment and subsequent reattachment of the protein 
at a point many base pairs away, and the basic quantity 
of interest is the time required for the protein to find its 
target on a finite DNA chain 0. Q The two mecha- 
nisms of diffusion and detachment/rcattachment can be 
viewed as a random walk along the chain with a variable 
step length distribution [j| ■ This is the viewpoint that 
we shall adopt for this work. A second motivation for our 
work is that when individual step lengths are drawn from 
a continuous distribution, the resulting random walk ex- 
hibit a variety of interesting properties beyond those of 
the discrete random walk. These include, among oth- 
ers, both unusual first-passage properties as well as 
extreme- value phenomena Q- 

Here we study the related problem of first-passage 
properties of the URW in a finite interval. Perhaps the 
most basic such feature is the exit probability i?(x), de- 
fined as the probability for a walk that starts at x to 
eventually cross the boundary at x = 1. Since the prob- 
ability to exit via the left boundary is 1 — i?(x), we need 
only consider exit to the right boundary, for which we re- 
serve the term exit probability. A related quantity is the 
mean exit time t(cc), defined as the average time to exit 
the interval at either boundary when the walk starts from 
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x. This exit time from an interval has been investigated 
for a random walk with general single-step hopping prob- 
abilities 0, E3 and also in an econophysics context 
0. As we shall see, the interplay between the maximum 
step size and the interval length leads to commensuration 
effects that are absent when the step length distribution 
extends over an infinite range. 

For both pure diffusion and the classical random walk 
with the step length Ax — l/N, where N is an arbi- 
trary integer, it is well known that the exit probability 
from an absorbing interval is R(x) = x Il2||. Simi- 
larly, the mean exit time is t(x) = x(l — x)/2D, where 
the diffusion coefficient for the discrete random walk is 
D = ((Ax) 2 )/2, with ((Ax) 2 ) the mean-square length 
of a single step. We now determine these first-passage 
properties — the exit probability and the exit time — for 
the URW. 



II. GENERAL FEATURES 

To help visualize the general behavior, we performed 
numerical simulations of the URW by a probability prop- 
agation algorithm; for the URW, this approach is orders 
of magnitude more efficient than direct Monte Carlo sim- 
ulation. In probability propagation, we first divide the 
unit interval into N discrete points. We correspondingly 
discretize the URW as follows: a URW that moves uni- 
formly within the range [—a, a] at each step is equivalent 
to its discretized counterpart hopping equiprobably to 
any one of the aN discrete points on either side of the 
current site. When an clement of probability hops out- 
side the interval, this element is considered to be trapped 
at the boundary where the element left the interval. We 
continue this propagation until less than 10 -6 of the ini- 
tial probability remains in the interval. Running the 
propagation further led to insignificant corrections. The 
simulations were generally performed with N = 5000. 
We found negligible differences in our results when the in- 
terval was discretized into N = 10000 points. This prob- 
ability propagation algorithm is also considerably more 
efficient than naive Monte Carlo simulation of an ensem- 
ble of random walkers. 
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When a is of the order of 1, non-diffusive features 
arise because the walk can traverse the interval in just 
a few steps. In the limit a — > oo, R(x) = 1/2 for all x; 
that is, either boundary is reached equiprobably, inde- 
pendent of the initial position. Conversely, for a <C 1, 
the first-passage properties of the URW approach those 
of continuum diffusion. Thus R(x) « x, except when the 
starting point is close to a boundary (Fig. In this 
boundary region, R(x) undergoes a series of transitions 
in which the n th derivative is discontinuous whenever x 
or 1 — x passes through na. These transitions become 
more apparent upon plotting R'{x) versus x (Fig. [21 . For 
small a, the qualitative behavior of R' is reminiscent of 
the Gibbs' overshoot phenomenon [Tflj when expanding 
a square wave in a Fourier series. 
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FIG. 1: Probability of exit at the right boundary, R(x), versus 
x for a = 0.9, 0.5, and 0.1 (slopes increasing, respectively), 
where a is the width of the single-step distribution. The data 
for a = 0.1 are based on probability propagation, while the 
other two datasets are obtained analytically. 

It is worth noting that the URW in an absorbing in- 
terval is not a martingale [lj| and hence R(x) ^ x in 
general. That is, the URW is not a "fair" process. The 
unfairness arises because when a walk crosses the bound- 
ary and ostensibly lands outside the interval, the walk is 
reassigned to be trapped at the edge of the interval. Thus 
the mean position of the probability distribution is not 
conserved. Another consequence of the unfairness is that 
the exit probability R(x — 0) for fixed a, no matter how 
small, is non-zero. The possibility of starting at x = 
and exiting at x = 1 can be viewed as an effective bias of 
the walk toward the middle of the interval. In the con- 
text of the gambler's ruin problem , a gambler that is 
about to be ruined appears to be best served by making 
a reckless bet that is of the order of the total amount of 
capital in the game. 

The mean exit time of the URW deviates strongly from 
the continuum diffusive form t[x) — x(l — x)/2D, when 
a becomes of the order of 1 (Fig.|3}. Notice also that t(x) 
does not go to zero as x — > or x — > 1. This limiting be- 




FIG. 2: The derivative R'(x) versus x for a = 0.9, 0.7, 0.5, 
0.4, 0.2, 0.1, and 0.05 (rising curves, respectively). Shown 
dashed are the results of an asymptotic approximation that 
is discussed in Sec. II VI 
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FIG. 3: Mean time to exit at either boundary, t(x), times a 2 
(proportional to the diffusion coefficient) versus x for a —* 0, 
a = 0.2, 0.4, 0.5, 0.7 and 0.9 (bottom to top). The data for 
a < 0.5 is based on simulation, while the remaining data is 
obtained analytically. 



havior again reflects the fact that there is a non-negligible 
probability for a walk that starts at one edge of the inter- 
val to exit via the opposite edge, a process that requires 
a non-zero time. 



III. FIRST-PASSAGE PROPERTIES 

In general, the exit probability may be determined 
from the backward equation that expresses the exit prob- 
ability R(x) in terms of the exit probability after one 
step of the random walk has elapsed 0] . This backward 
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equation has the generic form 

R(x) = / dx' p(x — > x')R(x') . 



(1) 



That is, the exit probability starting from x equals the 
probability of making a single step to x' times the prob- 
ability of exit from x', integrated over all possible values 
of x' . In a parallel fashion, the mean exit time can gener- 
ically be written as 



t(x) 



dx' p(x — > x )[t(x ) + 1] 



(2) 



That is, the exit time starting from x equals one plus the 
exit time from x', when integrated over all possible values 
of x' , with each term weighted by p(x — > x'). Note that 
the trailing factor of 1 can be taken outside the integral 
sum since J dx' p(x —> x') — 1. We now apply these 
two formulae to determine first-passage properties for the 
URW by studying, in turn, the cases a > 1, a G [1/2, 1], 
a G [1/3, 1/2], a < 1/3, and finally a -> 0. 



A. a > 1 

When a > 1, the support of the probability distribu- 
tion necessarily extends beyond the unit interval after 
one step. The residue that remains in the interval is also 
uniformly distributed. These facts allow us to obtain 
R(x) and t(x) by simple probabilistic reasoning. After 
one step, the walk jumps past x — 1 with probability 
x+ 2a~ 1 , an d jumps to the left of the origin with proba- 
bility (Fig. 0J . Because the remaining probability of 
l/2a is uniformly distributed within [0,1], there is a 50% 
chance that this residue will eventually exit via either 
end. Thus the exit probability is 



R(x) 



2a 



1 1 

2 2a~ 



2{x + a) - 1 
4a" 



(3) 



Since the first-passage probability for the walk to first 
exit the interval at the n th step is F(n) = S(n—1) — S(n), 
the mean exit time is 

oo 

(t) = £>[S(n-l)-S(n)] 
1 



n=l 

oo 

n=0 



1 - 1/2<Z ' 



(4) 



As a — > oo, (t) — * 1, while as a — > 1 from above, (t) — > 2. 
This same value for the exit time can also be obtained by 
solving the backward equation for t(x) itself (see below). 



B. a £ [1/2, 1] 

When 1/2 < a < 1, the unit interval naturally divides 
into an inner subinterval (1 — a, a) and outer subinter- 
vals (0, 1 — a), and (a, 1) (Fig. 0. If the walk begins in 
(1 — a, a), then the probability distribution of the walk 
necessarily extends beyond [0,1] after a single step, and 
the exit probability is again given by Eq. 
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FIG. 5: The unit interval for 1/2 < a < 1. A walk starting 
in (1 — a, a) can leave the interval via either end in a single 
step. 

On the other hand, when the walk starts in the outer 
subintervals the recursion formula Q for the exit prob- 
ability becomes 



h I dx'R(x') 
R(x) = { ° 

± J dx'R(x') 



x G [0, 1 - a] 
x G [a, 1] . 



(5) 



J \ 



x-a ox l x+a 

— a-x — — »!x+a-l^ — 

FIG. 4: Probability distribution of the uniform random walk 
after a single step. The shaded regions are the portion of the 
probability distribution outside the unit interval. 

As expected, the exit probability approaches 1/2, inde- 
pendent of the starting point, as the average step length 
becomes large. Also notice that R(l) = \ + asa^ 1 
from above, the probability of exiting the right boundary 
when starting at x — 1 is only 3/4. 

The survival probability for the walk to remain within 
the interval after a single step is simply l/2a, and the 
survival probability after n steps is then S(n) = (l/2et)™. 



Differentiating these equations gives 



2aR'(x) = R(x + a) 
2aR'(x) = 1 - R(x - a) 



x G (0, 1 - a) , 
x G (a, 1) . 



When x G (0, 1 — a), the second derivative is 
4a 2 R" (x) = 2aR'(x + a) = 1 - R[x), 



(6a) 
(6b) 



(7) 



where we use the fact that if x G (0, 1 — a), then x + a E 
(a, I). The solution to J7J is 



R(x) = 1 + Ci sin — (x — C2) 
.2a 



(8) 



To determine the constants c\ and C2, we first substi- 
tute JSJ into (|6a|) and also use the fact that R(x + a) = 
l — R(x) for x G (0, l-a) to obtain c 2 = (l-a)/2 + 7ra/2. 
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Second, we match (|HJ and Q at x — 1 — a to find C\. 
The exit probability for x £ (0, 1 — a) therefore is 



R(x) 



1 + 



1 

2a 



1 -a 



a- 1 
4a 



(9) 



The sinusoidal segment of R(x) is visually close to a lin- 
ear function (Fig. Q), and the difference between these 
two functional forms becomes more clearly visible upon 
plotting R'(x) versus x (Fig. ^J. 

We now compute the mean exit time. Again, there 
are two cases to consider: either the walk begins within 
(a, 1 — a) or it begins in the complementary outer subin- 
tervals. Let us denote by t- lVL (x) and t ou t(x) as the mean 
exit times when the walk starts at a point x in the in- 
ner and in the outer subintervals, respectively. Then the 
backward equation @ for t- m becomes 



tin (a;) = 



1 

1 + — / t(x')dx' 
2a J 



o 

1 — a a 

1 + ^ f t aat (x')dx' + ^ Jt iu (x')dx'.(lQ) 

1-a 



For the last line, we break up the integral into a contribu- 
tion from the outer subinterval, with two equal contribu- 
tion from (0, 1 — a) and (a, 1), and the inner subinterval 
(1 — a, a). Notice also from the first line that ti n (x) is 
independent of x. Thus we define t- ln (x) = T, with T 
dependent only on a. 

Similarly, the backward equation for t out is 



*out(a;) = 1 + 



1 

27i 



x-\-a 



t(x') dx' . 



(11) 



Differentiating gives t out (x) — t out (x+a)/2a. Notice that 
if x £ (0, 1 — a), then x + a is necessarily in (a, 1). Cor- 
respondingly, the backward equation for t out (x + a) gives 
t' out (x + a) = -t out (x)/2a. Thus % ut (x) = -t out (x) / Aa 2 , 
with solution 



tout (a;) = n cos ( — ) 
V la) 



r2Sin (£ 



(12) 



To complete the solution, we need To determine the 
three unknown constants T, n, and r%. The solution is 
straightforward and the details are given in Appendix 
From this solution, quoted in Eq. (|A6|> . we obtain the 
mean exit times plotted in Fig. |3] As a decreases, t(x) 
quickly approaches the parabolic form of the diffusive 
limit, but t(0) and t(l) = t(0) remain strictly greater 
than zero when a is non zero. 



C. a e [1/3, 1/2] 

For any a < 1/2, the exit probability now obeys the 
generic recursion formulae: 



i Tdx'R(x') 
o 

= { k Tdx'R(x') 

x — a 

± J dx>R(x>) 



,r • (i — 1 



x e [0, a] 

x € [a, 1 — a] 
x G [1 - a, 1]. 



(13) 

For example, the middle equation states that the exit 
probability starting at x equals the exit probability after 
making one step to x' — which is uniformly distributed 
in the range — [a, a] about x — times the exit probability 
from x' . The first and third equations account for the 
modified range of the single-step distribution if the walk 
leaves the interval. 

We differentiate these integral equations to obtain the 
more compact form 



2aR'(x) = R{x + a) 
2aR'(x) = R(x + a)~ R(x 
2aR'{x) = 1 - R(x - a) 



x e (0,a) ; (14a) 
x G (a, 1 — a) ; (14b) 
x G (1 - a, 1) . (14c) 



To solve these equations for the cases where o G 
[1/3, 1/2], we should consider the five subintervals (0, 1 — 
2a), (1 - 2a, a), (a, 1-a), (1-a, 2a), and (2a, 1) (Fig.®. 
By symmetry, we only need to study the range x < 1/2, 
and we now examine, in turn, subintervals III, I, and II. 
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FIG. 6: The unit interval with the subregions used to deter- 
mine the exit probabilities for the case a £ (1/3, 1/2). 



Subinterval III: For x £ (a, 1 — a), Eq. (|14b|l con- 
nects R' (x) with R(x + a) and R(x — a). In turn, the 
equations for R'(x + a) and R'(x — a) involve R(x) with 
x again within (a, 1 — a). Thus, 

Aa 2 R"(x) = 2aR'(x + a)-2aR / (x-a) = l-2R(x) , (15) 

with solution 

1 ■ T 1 ( 1\" 
— h Ci sm — =— x 

2 [V2a\ 21 

This form automatically satisfies the symmetry condition 
J2(l/2) = 1/2. 

Subinterval I: We obtain the exit probability for x £ 
(0,1 — 2a) by integrating Ijl4a(l and also using the fact 
that the argument x + a in R(x + a) lies within (a, 1 — a) . 
Thus we use the result of Eq. (fTB|> to give 



R(x) 



x £ (a,l 



(16) 



R(x) = = cos 

4a V2 



1 



V2a 



x+a- 1 - 



-c 2 x £ (0, l-2a) 
(17) 
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To determine the constants c\ and C2 , we use the general 
antisymmetry condition, R(y) = 1 — R(l — y), to write 
(I14bll in the form 



near x = 2a. By l|14b|) . i?"(ir) is coupled to R'(x + a) 
and R' (x — a), and the latter derivative has a jump when 
its argument equals a. Thus R"(x) has a jump of magni- 
tude i?(0)/4a 2 when x passes through 2a. This pattern 



2aR'(x) = l-R(l-x~a)-R(x-a) x G (a, 1-a) . (18) continues so that for a E ~), the n th derivative of 

Now we substitute the solutions l|16(l and H17JI into 
Eq. O and find c 2 = | - §k Thus 



1/2 ci 
V2 



4a 



cos 



/2a V 



\/2a V 2 



(19) 



Subinterval II: Finally, for x G (1— 2a, a), Eqs. Ijl4a(l 
and H14cfl show that the subintervals (1 — 2a, a) and 
(1 — a, 2a) are coupled only to each other. Using the 
antisymmetry of the exit probability about x = 1/2, we 
have 



R(x) has a jump discontinuity as x passes through na, 
while all lower derivatives are continuous. Thus R(x) be- 
comes progressively smoother and more linear in visual 
appearance for x deeper in the interior of the interval. 
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FIG. 7: The unit interval (to scale) for the case | < a < j 
and the two associated classes of subintervals. The shaded 
subintervals are labeled sequentially as defined in the text. 



4a 2 R"(x) = 1 - R(x) x€(l- 2a, a) 



(20) 



with solution R(x) = 1 + C3 sin((a; — C4)/2a). We deter- 
mine C4 from the condition R(x) = 1 — R(l — x) to then 
give 



R(x) — 1 + C3 sin 



1_/ 1-a 
2a~r ~^2~ 



(21) 



To obtain the remaining two constants ci and C3, we 
match Eqs. JE)) and J2T1) at 1 - 2a, and J2T1) and JTHJl at 

a. These lead to 



ci 



C3 



• tana 



-4= cos — sin 8 tan a 

V2 



J_ 

8a 



■ tan a 



2 cos a cos a cot /3 — sin a ' 



with 



3a 



4a 



7T 

4' 



1 - 2a 
2V2a ' 



(22a) 
(22b) 

(23) 



For the special case of a = 1/3, subinterval II disappears 
so that the solution consists of (|17|) and l|16fl only, and 
the relevant constant in l|22a|l simplifies to 



(24) 



c\ = [ 4\/2 cos —= — 8 sin —= 



D. 



a < 1/3 



It is straightforward to treat smaller values of a, but 
the bookkeeping of the various subintervals becomes in- 
creasingly tedious. However, it is still possible to in- 
fer general properties of the exit probability. From 
Eqs. (|Hajl & ifilb) . we have 2aR'(x = a + ) = R(2a) - 
R(0) while 2aR'(x = or) = R(2a). Thus R'(x) has a 
jump of magnitude i?(0)/2a when x passes through a, as 
illustrated in Fig. [21 a < 0.4. Similarly, consider R"(x) 



The analytical solution for R(x) can, in principle, be 
obtained from the backward equations (|14afl - l|14cf> . for 
the exit probability. When a < 1/2, these equations 
naturally partition the unit interval into two classes of 
subintervals as shown in Fig. [3 Nearest-neighbor shaded 
subintervals are coupled only to each other by these back- 
ward equations, and similarly for the complementary 
subintervals. It is convenient to define r^{x) as the exit 
probability for a walk that begins at x + (k — l)a in the 
k th shaded subinterval; that is r^ipo) = R(x+(k—l)a) for 
x £ {(k— l)a, ka). With these conventions, the backward 
equations for the exit probability for a starting point in 
one of the shaded or in one of the unshaded subintervals 
have the matrix form: 



2ar'(a;) 



/ 1 

-1 1 
0-10 



V 



0-1 01 
0-10/ 



r(x) 



(25) 



where v(x) is the column vector with components 
[r 1 {x),r 2 {x) . . . r n (x)} = [R(x),R(x+a) . . . R(x+(n-l)a)} 
and I is the column vector with the n components 
(0,0,. ..,0,1). 

To solve these equations, we note that the eigenvalues 
of the above matrix are given by Aj = 2i cos(nj/ (n + 1)), 
with j = 1,2,..., n Thus r n (x) is, in general, a 

linear superposition of the eigenvectors of the matrix; 
these are sinusoidal functions with arguments Xj/2a. The 
actual form of R(x) is then obtained by fixing the various 
constants in this eigenvector expansion through matching 
the components at appropriate boundary points. 



IV. EXIT PROBABILITY FOR a 







In the limit a — > 0, the diffusion approximation be- 
comes increasingly accurate so that R(x) is very nearly 
equal to x, except in a small region of the order of a near 
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each boundary (Fig. ^) . This deviation is more clearly 
evident when plotting R'(x) versus x (Fig. As a gets 
small, this plot also suggests that a good approximation 
to R'(x) will be obtained by solving the exact equation 
for R'(x) in the boundary region and treating R'(x) as a 
constant in the interior of the interval. 

At a zero-order level of approximation, we assume that 
R'(x) is position dependent for x S (0, a) and x G (1 — 
a, 1) and is constant otherwise. Then within (0, a), the 
backward equation 2aR'(x) — R(x + a) means that R'(x) 
is a linear function (and similarly for x G (1 — a, 1)). Thus 
we make the ansatz 



R{x) 



|(1 — s) + sx x G (a, 1 — a) 
tq + r\x + T2X 2 otherwise , 



(26) 



with s and the ri to be determined. We expect the slope 
s of R(x) in the interior of the interval to approach 1 as 
a — > to recover R(x) — > x in this limit. The form of the 
constants in the first line also ensure the obvious special 
case R(l/2) — 1/2. Similarly, the linear form for R'(x) 
in the boundary regions roughly corresponds to what is 
seen in Fig. [3 




FIG. 8: Exit probability R(x = 0) versus a, based on 5 x 10 8 
walks for 200 equally-spaced values of a between and 2. 
Shown dotted is the result of our zeroth-order asymptotic 
approximation -R(O) = a/[2(l + a)], while the dashed curve is 
the result of the first-order approximation. 

We determine the 4 unknowns in the above asymptotic 
approximation for R(x) by the following conditions: (i) 
2aR'(x) = R(x + a) must be satisfied in the region x G 
(0, a) (this gives two conditions — one for the linear term 
and one for the constant term), (ii) R(x) is continuous, 
(iii) the discontinuity in R'(x) at x = a equals R(0)/2a, 
as follows from Eqs. i|14a|) and i|14b(l . Applying these 
conditions gives, after some simple calculation, 



1 



1 



3 1 



1 



2 1 



r l = T 



4 1 



T2 = 



1 1 

4a~TT 



A much better approximation is obtained by treating 
R(x) exactly in the domains (0, a) and (a, 2a), and then 



assuming that R'(x) is constant otherwise. Thus for x > 
2a, R(x) is given by the first line of Eq. (|26[) . while within 
(0, a) and (a, 2a), the governing backward equations for 
R(x) are 



2aR'{x) 



R(x + a) x < a 

R(x + a) — R(x — a) x > a 



(27) 



For x < a, we iterate the first equation to give 
Aa 2 R"(x) = R(x + 2a) — R(x) and make use of the as- 
sumption that R'(x) = s for x > 2a. This leads to the 
approximation 

d 2 sin (|-) + s(a;+2a) + = (1-*) . 

(28) 

Similarly, to obtain R(x) in the region (a, 2a), we in- 
tegrate the backward equation 2aR'(x) = R(x + a) — 
R(x — a) and again use the fact that the argument x + a 
in R(x + a) is beyond 2a, so that R(x + a) is a linear 
function. This integration leads to 



R(x) — di cos ( — 
\2a 



R(x) = —d\ sin 



2n 



G?2 COS 



2o 



(29) 



We determine the 4 coefficients in these two forms for 
R{x) by requiring that at x — a, R is continuous and 
the first derivative has a jump of magnitude i?(0)/2a, 
while at x = 2a, both R and R' are continuous. The 
resulting formulae are given in Appendix iBl Fig. [21 shows 
the result of this small-a approximation for R'(x). The 
agreement between this asymptotic approximation and 
the numerical results is extremely good. 

As a further test of the accuracy of this approach, we 
show the numerically-obtained dependence of R(0) on a 
together with our zeroth- and first-order approximations 
for R(0) (Fig. EJ|. As already mentioned in Sec.|TT| R(0) 
is greater than zero for any a > because there is a non- 
zero chance that a walk exactly at the left boundary can 
still exit via the right boundary. In the limit of small a, 
our asymptotic approach for i?(0) closely approximates 
the data. 



V. CONCLUDING REMARKS 

An enigmatic feature of the uniform random walk 
(URW) is that its first-passage properties in a finite in- 
terval are not described in terms of a radiation boundary 
condition [l2L ITrjl ]. This boundary condition is d = c/k, 
where n is the extrapolation length. This condition can 
be interpreted as partial absorption and partial reflec- 
tion when a walk hits the boundary. In the case of a 
URW that starts at one absorbing boundary, there is 
a non-zero chance for the walk to exit via the opposite 
boundary. This incomplete absorption should be equiva- 
lent to partial reflection from the initial boundary, which 
should be described by a radiation boundary condition. 
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Indeed, the probability distribution of a URW in a semi- 
infinite interval x > with absorption at x — closely 
approximates that obtained for pure diffusion with a ra- 
diation boundary condition. However, in the finite inter- 
val the radiation boundary condition gives R{x) as the 
linear function R(x) = y+Tk which does not account for 
the anomalous behavior observed near the edges of the 
interval. 

What we do find is that the first-passage properties 
of the URW in a finite interval exhibit curious com- 
mensuration effects as a passes through 1, 5, g ■ ■ ■• For 
a e (^ji ~), the n th derivative of R(x) has a jump dis- 
continuity as x passes through na, while all lower deriva- 
tives are continuous. The exit time has corresponding 
singular behaviors. For small n, we have computed the 
exit probability and the mean exit time exactly by a 
direct probabilistic approach. In the limit a — > 0, the 
diffusion approximation becomes increasingly accurate, 
except when the starting point is close to either bound- 
ary where the exit probability continues to exhibit non- 
diffusive effects. In the limit of small a, we constructed 
an approximation of treating the exit probability exactly 
within the boundary region, an approach that gives ex- 
tremely accurate results. 
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For Eq. Qllfl- we evaluate it at x = 0. This gives 
I 



tout(0) = I + 77- / t(x')dx' . 
2a Jo 

Then using Eq. (fT2"|) for f O ut(0) and separating the inte- 
gral into the inner and outer subintervals, we have 



71 = 1 + S T+ ( 1 -S' :r - 



(A4) 



Finally, we equate t m and t out at x = 1 — a, where the 
inner and outer subintervals meet. This gives 



7 = Ti cos 



1 - a 
2a 



t 2 sin 



1 - a 
2a 



(A5) 



The conditions l|A2|l . (|A4(I . and l|A5(l provide the three 
independent equations 



1 1 — a , 

7 = 2a + 4an sin — — - 4aT 2 



2a 



ti = 1 + an sin — — 1 - r 2 



1-a , 
cos I — — — I 



2 a 



2 a 



I — a 
cos I — — — 1 



2a 



1— - 



2 a 



T = n cos 



1-a 



T 2 COS 



1-a 



2a J \ 2a 

for the unknown coefficients T, n, and t 2 . 

To express the solution succinctly, let z = a 
e = (l — j-). Further, define 

a — 4a sin z — cos z [3 — 4a(cosz — 1) + sin 2; 
7 = sin z — 1 + e cos z 8 = cos z — 1 — e sin z . 

Then the constants are 



APPENDIX A: MEAN EXIT TIME FOR a > 1/2 

To complete the solution for the mean exit time, we 
substitute t m +a and i out quoted in Eq. $T2$ into Eqs. ijTU|) 
and l|llfl . The former equation becomes 



T= 1 



1 



i ou t(z')aV + ^r^ T - ( A1 ) 

2a 



This can be rewritten as 

7= 2a + 2T, 

where 

pi — a 

T = t out (x')dx' 
Jo 



2a 



2 a 



(A2) 



TiC0S {k) +T2Sin {h) dx ' 

= 2an sin ( — - — ) — 2ar 2 cos ( — - — ) . (A3) 



2aS-P 2a + an 
n = — r 2 = 7 = Ti COSZ + T2 smz.(A6) 

-ffj — ao p 

The result of this solution is shown in Fig. [3] 



APPENDIX B: COEFFICIENTS OF THE EXIT 
PROBABILITY FOR a -> 

At x = a, R{x) is continuous, while R'(x) has a dis- 
continuity of size R(0). Similarly, at x = 2a, both 
and R'(x) are continuous. From Eqs. (|28|l and l|29|l . we 
thus have the conditions: 

di cos i + d 2 sin - + 3.sa + -(1 — s) = d 2 + d3 

1 11,, 
-di sin - + d 2 cos - = - (1 - s) 

—d\ sin 2+^2 cos — + d% = +2sa + — (1 — s) 

— g?i cos — — d 2 sin — = 2sa . (Bl) 
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The solution to these equations are: with 

(l-«)/2 



di 



d 2 



( u - ^(s 008 ! + sin ^) + ( cos I - s sin §)« 

di 1 d 2 . 1 111111,, 
— — cos - — — sin - u = cos - + - sin - v = — cos - — sin - . (B3) 

2a 2 2a 2 222222 v/ 

~ u- 1 



1 11 

d 3 = 2sa + -(1 - s) + d\ sin - — d 2 cos - , (B2) 
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